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Abstract 

We study the effect of hole hopping in a doped antiferromagnet described by 
the 'i-J ' model, using exact analytic solutions for small clusters. In spite of 
the small size, they reveal interesting details about the magnetic order, which 
are not apparent in Mean Field treatments or in numerical calculations. The 
4-site cluster with one hole yields the most interesting physics, displaying 
different behaviors for the ground state: an antiferromagnetic phase for 
t <C J, where the hole seems to be localized, not affecting the order of the 
Heisenberg spins; ii) another regime for t ^ J that presents mixed ferro and 
antiferro correlations and coexistence of metallic and insulating behaviors, 
with the presence of charge and spin density waves, in what may be the 
analog of the spiral phase obtained in Mean Field solutions; and finally Hi) 
for t ^ J , we obtain strong ferromagnetic correlations (maximum spin) and 
no density waves, with quantum fluctuations precluding the saturation of the 
magnetic moment. This behavior shows traces of a metal-insulator transition 
as the hole kinetic energy competes with the antiferromagnetic interactions. 
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The physics of the normal state of High-Tc compounds is yet to be understood after nearly 
12 years of great experimental and theoretical developments. The antiferromagnetic long 
range order at T = is rapidly destroyed when the CuO planes are slightly doped with holes 
(x ~ 0.02 for La2-xSrxCu04) [Q]. These holes form 0~ ions which coupled to the Cv?^ spin 
1/2 local moment can be regarded as singlets centered at the coppers. Looking at this singlet 
as a hole on a Cu^'^ square lattice, Zhang and Rice derived an effective Hamiltonian, 
the so called 't — J ' model, which describes the hopping of holes on a spin lattice, with the 
spins coupled via an antiferromagnetic Heisenberg interaction. Mean field theory (MFT) 
was extensively used by numerous authors to treat the hole motion on the interacting spin 
background 0-^. They found an evolution from an antiferro to a ferromagnetic phase, 
when the parameter t/J is increased at low or intermediate hole concentrations. The most 
interesting feature of those calculations is the spiral phase that they predict for t ~ J, a trait 
that seems to have been observed in La2-xSrxCu04 by neutron scattering experiments 0. 
The question then arises on what would be the effect of quantum fluctuations (which are 
neglected in MFT) on this magnetic order, and what is the nature of this spiral state. Since 
this phase is induced by hole mobility, it should be accompanied by an insulator-to-metal 
transition. 

In this paper we calculate exact eigenf unctions for the 't — J ' model Hamiltonian 

Ht-j = -t 5^ Po (clc,. + h.c.) Po + JJ2Po ( 5^. ■ 5^. - ^) ^o, (1) 

cr<i,j> <i,j> 

where is the fermionic operator that creates an electron (hole) with spin a on site i, 
rii = ^la^io- is the number operator, Pq projects out double occupied states in the lattice, 
and the brackets (...) under the summation signs mean nearest neighbors. The parameter 
t measures the mobility of carriers which hop from site to site in an antiferromagnetic 
Heisenberg background, with exchange J > 0. We were able to obtain exact analytic 
solutions for the eigenf unctions of (|^) in a 4-site cluster, which despite its small size, shows 
traces of a metal- insulator transition for the ground state (T = 0), when the parameter 
A = t/ J is varied. We also propose a way to analyze the subjacent magnetic order for A ~ 1, 
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searching for a precursory spiral phase in this small cluster. We solve the problems with one 
and two holes on a square cluster. The two-hole case is simpler, since the remaining spins 
accommodate to form a Resonating- Valence-Bonding state for small A, and holes develop 
Charge Density Waves (CDW) states. The other phase (A 1) is a spin triplet and the two 
holes are delocalized. The physics of the one-hole problem is much richer, since there is a 
spin that cannot be compensated in spite of the antiferromagnetic correlations. As a result, 
we obtained an intermediate phase, with mixed insulating and metallic behavior, displaying 
a variety of magnetic textures, one of which we identify with the spiral phase predicted by 
MFT. We are aware that a small cluster cannot display true phase transitions. However, 
broken spin and charge symmetry states can be constructed for the ground level that, at 
least qualitatively, can be associated with macroscopic phases. [0] 

We first consider the case of two holes on a square cluster. This problem has dimension 
6x2^ = 24. The ground state energy shows a crossover at A = ^^^ZT) - The small A 

phase has an antiferromagnetic character, with vanishing total spin and total component 
(note that commutes with and j for periodic boundary conditions). It is 

double degenerate and can be written as a dimerization of singlets on the cluster, as shown 
in Fig. 1, where the double links stand for spin singlets and the lone circles for holes. This 
resembles closely the Resonating- Valence- Bond state (RVB), as prescribed by Anderson for 
the Heisenberg model Broken symmetry states in the form of Charge Density Waves 
(CDW) can be obtained by linear combinations of the above (for example, (ki > +|cf2 >) 
localizes the holes on one of the diagonals). The above features are characteristic of the 
insulating state. In contrast, the ground state for A > 1.21 shows a qualitatively different 
behavior: it is a spin triplet {Sz = 1, 0, —1, see the states below), with admixtures from 
all the 24 kets of our base, suggesting some type of disorder induced by hole hopping. They 
are given by: 
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with an obvious notation for the basis functions that depicts the geometry of the square, 
the symbol standing for holes. Note that basis kets entering in one of the states \bi) do 
not appear in any other meaning that linear combinations of these degenerate states 
will not lead to CDW's (the probability of measuring a hole in one of the sites is the same 
for all sites). Therefore the states can be seen as precursors of the metallic regime. It 
comes as a surprise that we can see traces of a metal-insulator transition on such a small 
cluster. 

Now, we turn to the problem of four sites with one hole, which is closer to the half 
filled band condition and displays more interesting magnetic properties. Before obtaining 
the exact ground state, we propose a variational solution for the case of low hole mobility 
(A ~ 0). This way, we expect to shed light on approximations for larger cluster problems, 
since analytic solutions can only be obtained in very special cases. As shown below, it turns 
out that our variational state corresponds to the exact solution when A is smaller than a 
critical value. For low t (or high magnetic coupling J), we expect the magnetic order to 
resemble that of the Heisenberg model for a 3-site chain with open boundary conditions. 
The ground state for this latter problem has energy — J and is 

1 



i^?„) = — {iTn)-2inT) + iiTT)}. 



(2) 



Next, we define the state \(f)j) as having a hole on the j-th site with the magnetic order given 
by (Of) for the remaining sites, and propose the variational wave function 
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i=i 

suggesting that the hole motion just introduces a phase 6j as it resonates through the cluster. 
We find the mean energy E = (^| Ht-j |\E') to be 

E = -J+^ [cos (^1 - 62) + cos {02 - 93) + cos (^3 - Oi) + cos (^1 - ^4)] , 

which is minimal for 61 = , 62 = tt , 63 = 27t , 64 = Svr, with the value — J — t. For what 
range of A is this variational solution valid? We will show below that for < A < |, this is 
indeed the exact ground state for the cluster. This state is degenerate with the one obtained 
by reversing all the spins, in the negative sector of the magnetization. 

To get the exact solution, we note that the Hilbert space has dimension 4 x 2^ = 32 (the 
hole can occupy each of the 4 sites and each spin can be up or down), but the Ht_j matrix 
divides up in 2 blocks of 5"^ = ±| (each of dimension (4x4)) and 2 of 5*2 = ±| (each 
of dimension (12 x 12)). We only diagonalize the positive magnetization blocks, the other 
eigenfunctions being obtained by fiipping the spins of the basis (from now on, we always 
work on the positive magnetization sector). In addition, the (12 x 12) matrix, when written 
in the total spin basis, breaks up in 2 blocks of 8 x 8 (spin 1/2) and 4x4 (spin 3/2). The 
energy levels can be seen in Fig. 2 and the eigenstates are given in the Appendix. The 
first feature we note in Fig. 2 is the great number of level crossings at J = 2t. This is the 
Supersymmetric point, where an exact solution via the Bethe Ansatz has been obtained in 
one dimension [^, with exact separation of the charge and spin degrees of freedom [jlO 



We see that the ground state energy undergoes two level crossovers, with the total spin 
being 1/2, 1/2 and 3/2, as the parameter A increases. This cascade effect on the total spin 



has already been reported in numerical calculations |TT],|T2||. The first ground state function 
holds for < A < i, has energy E = —J — t, and it is exactly equal to our variational 
solution shown above. The hole motion does not affect the magnetic order, just changes the 
relative phase between antiferromagnetic components. Although we expected CDW's, this 
expectation is not realized, since there is not enough degeneracy in the ground manifold. 



The absence of a broken symmetry state might be ascribed to the small size of the system. 
In this case, size effects are lifting the degeneracy with the first excited level, degeneracy 
that is restored at the Supersymmetric point. 

The second ground state is double degenerate (|ci) and |c2) in the Appendix) and holds 
for I < A < ^ 3.8, with energy E ^ -^(1 + ^/lTW^). What makes those states 

non trivial is the fact that they depend on A (see the functions a,b, c, ... in the Appendix), 
making them harder to analyze. Both of them obey a curious rule of sign when the hole 
hops, the same pattern repeating with a negative sign, in the same way as in the variational 
solution. It is easy to see that we can build linear combinations of |ci) and |c2) (which are 
linearly independent but not orthogonal), with real coefficients, that concentrate the hole 
on the diagonals, i.e. CDW's can be obtained with those states. Fig 3 shows the probabihty 
Pi of finding the hole at the i-th site as A varies, for the two kets |ci) and |c2). In the high 
mobility limit (A — > oo), all pi 's go asymptotically to 1/4, as expected. 

As noted in the introduction, MFT's yield a spiral magnetic order for A ~ 1. We would 
like to test if our solution resembles this spiral in any way. We note that the states |ci) and 
|c2) have no uniform spin distribution. Broken spin symmetry states in the form of Spin 
Density Waves (SDW) can be build with the above states in the form 



where ± refer to the manifolds Sz — ±1/2, respectively, and {9, (f) are arbitrary parameters 
which can be related with the geometry and pitch of the density wave. As noted before, 
states |ci±) and |c2±) are degenerate and linear independent and were chosen for their 



spin distribution within the cluster, with a net ferromagnetic component that precesses with 
0, and whose amplitude is modulated as A and 9 vary. The detailed calculation yields: 



\SDW) = cos 




ci+) + exp (i0) sin 




simple pattern. Calculation of the site magnetization (SDW ~§ j SDW) shows a spiral 



< ^1 >=< 5^3 >= [h{X) cos 9 - z{X)]z- g{X) sm9 n± 

< >=< ^4 >= [j{X) cos 9 + z{\)]z + k{\) sm9 n± 



(3) 
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where n± = cos((/))x + sin(0)y is a unit vector on the xy-plane and {h, z, g,j, k) are functions 
given in the Appendix. Fig. 4 depicts this distribution for various choices of (0, 0). Some of 
the configurations shown closely resemble the double-spiral state found in Ref. through 
a mean-field calculation. Canting of the spins in doped one-dimensional antiferromagnets 
has also been reported previously by one of the authors using a variational approach []T3 



The presence of CDW and SDW in the above states suggests that they are precursors of the 
insulating regime. What is a remarkable peculiarity of the region | < A < is that 

we can also construct a ground state with uniform charge and spin distributions, i.e. with 
the same probability pi of finding the hole in any site and correspondingly, with the same 
site magnetization. We call them Uniform Hole Probability (UHP) states. Interestingly 
enough, the only two combinations that yield this result are complex and are only feasible 
for A > 1/2 : 

\uhp a) = Vl-u^ Ici a) ± iu Ics a ) , with u(X) = 12A2 - 1 - Vl + 12X^ , (4) 

6A ^ 



where the |...) symbol stands to indicate normalized states and a = ± refers to the manifolds 
^2 = ±1/2, respectively. For these states, we get 



Pi = P2 = P3 = P4 = i , 



(^.(l)) = (5.(2)) = (5.(3)) = (5.(4)) = a|. 
We can generalize the above state to 

\uhp n) = cos ^-^ \uhp+) + exp sin ^-^ \uhp—) , (5) 

as a state with uniform charge distribution and with the uniform magnetization pointing 
along the direction n (n being a unitary vector with latitude 6 and azimuthal angle 0). 
Clearly, the \uhp) kets are precursors of the metallic state. Thus, in the region | < A < 
beginning at the Supersymmetric point, the system displays a coexistence of metallic 
and insulating behaviors. Two main questions remain: i) why is the coexistence region so 
wide; and ii) how can we extrapolate these results to the infinite-size limit? The degeneracy 



of the ground state can be ascribed to the vicinity of the Supersymmetric point (J = 2t) and 
to the square symmetry (in a sense, it may be a particular feature of the 4-site problem). The 
large region spanned by this precursory phase is certainly a size effect, and it will extrapolate 
to a single point for the infinite-size limit, or to a line when the doping is introduced as a 
variable. 

The third ground state, for A > 3.8, presents the same metallic character as in the two- 
hole problem solved at the beginning of this paper: it presents ferromagnetic spin correlations 
but holes are delocalized, all the base kets having the same admixture, thus avoiding the 
saturation of the magnetic moment (this is clearly due to quantum fluctuations, contrary 
to what is observed in MFT's). This state is double degenerate and its energy is ^ — 2t. 
CDW's cannot be obtained, as in the two-hole case. 

In conclusion, we described traces of a metal insulator transition in the ground state as 
we increase A = j. Although this treatment can only be done in small clusters, we were 
able to do a clear analysis of the wave function, showing properties that are not apparent in 
a number of numerical calculations performed for the same problem |]r2|. We analyzed the 
evolution of the magnetic properties by explicitly writing the wave functions, in contrast to 
the usual self-consistent calculation of magnetic averages. Broken charge and spin symmetry 
states are proposed as precursors of macroscopic phases. In particular, we found Spin Density 
Wave states that can be identified with the spiral phase obtained by Mean-Field procedures. 
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A. Appendix: Eigenfunctions for the 4-site cluster with 1 hole 

The basis for — ^ is (note the order): 
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FIGURES 

FIG. 1. Bi-degenerate ground state for t/JK 1.21. The double links are singlets which 
characterize this state as a RVB insulator. 

FIG. 2. Energy levels of the 4-site cluster with 1 hole. As we increase A the ground state energy 
undergoes two level crossings. Note the amount of degeneracy for A = 0.5, the Supersymmetric 
point. 

FIG. 3. Probability of measuring the hole on the \-th site; solid line is for \c\) and dotted for 
|c2). Those states span the ground manifold when 1/2<A<3.8. 

FIG. 4. Spiral magnetic order of SDW states for t ^ J. In all cases, the angle ^ gives the 
orientation on the xy-plane: a) ^ = 7r/2; h) 9 = tt; c) 6 = 37r/2. The thin ends of the spins are 
meant to point into the page. 
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